Two different kinds of rogue waves in weakly-crossing sea states 
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Formation of giant waves in sea states with two spectral maxima, centered at close wave vectors 
ko it Ak/2 in the Fourier plane, is numerically simulated using the fully nonlinear model for long- 
crested water waves [V. P. Ruban, Phys. Rev. E 71, 055303(R) (2005)]. Depending on an angle 
between the vectors ko and Ak, which determines a typical orientation of interference stripes in the 
physical plane, rogue waves arise having different spatial structure. HO < arctan(l/V2), then typical 
giant waves are relatively long fragments of essentially two-dimensional (2D) ridges, separated by 
wide valleys and consisting of alternating oblique crests and troughs. At nearly perpendicular ko 
and Ak, the interference minima develop to coherent structures similar to the dark solitons of the 
nonlinear Shroedinger equation, and a 2D freak wave looks much as a piece of a ID freak wave, 
bounded in the transversal direction by two such dark solitons. 

PACS numbers: 47.35.Bb, 92.10.-c, 02.60.Cb 



The problem of extreme ocean waves (known as freak 
waves, rogue waves, or killer waves) has attracted much 
attention in recent years (see, e.g., the reviews [l|, 
where different physical mechanisms of the rogue wave 
phenomenon are discussed, and many related works are 
referenced; for some recent developments in this field, 
see Refs. 0, i, i, H, 0, i, B 113, El, E, El, Q E, 
fl6l . 1171]). With a typical background wave amplitude 
Aq ^ [0.015 ••• 0.02] Ao (where Aq = 27r//co is a typical 
wave length) , the maximum elevation of a freak wave can 
reach a height Fmax > O.O6A0, which approaches the lim- 
iting Stokes wave. Profiles of freak waves are very steep, 
and they strongly deviate from the sinusoidal shape. In 
different circumstances, the giant waves can be caused 
by different reasons. Accordingly, there are several prob- 
able scenarios explaining formation of these waves. It 
has been recognized that one of the most important rea- 
sons for freak waves is the modulational Benjamin-Feir- 
Zakharov instability taking place in relatively long and 
high groups of propagating waves Efficiency of 

this mechanism is usually characterized by the so-called 
Benjamin- Feir index (BFI) [4], 



BFI ^ Ao Mo/o, 



(1) 



where Iq is a typical length of wave groups. For exam- 
ple, in completely incoherent sea states (low BFI) the 
modulational instability is suppressed, and rather rare 
appearance of anomalous wave events is basically of a 
purely kinematic origin. This limit is well described 
by the approximation of non-interacting normal wave 
modes, renormalized by a weakly- nonlinear transforma- 
tion excludinff^ three- wave (non-resonant) processes (see, 
e.g., Refs.jil, m, ll^, and references therein). Higher val- 
ues of BFI correspond to more favorable conditions for 
the occurrence of freak waves. Nonlinear wave interac- 
tions become essential, so giant waves arise in the process 
of evolution of some coherent structures. In particular, 
the limit of infinitely high BFI (a weakly disturbed pla- 
nar wave as an initial state) has been recently studied 
in works [1,0], and specific zigzag-shaped, obliquely ori- 



ented wave stripes were found to develop in the nonlinear 
stage of the modulational instability, with rogue waves 
occurring mainly at zigzag turns. That case roughly cor- 
responds to another probable scenario, when refraction of 
swell in a spatially non-uniform current causes significant 
preliminary amplification of wave height around caustic 
region [2l|, [22[ . A different kind of coherent structures 
has been recognized recently for purely one-dimensional 
(ID ) waves (planar flows), the so-called giant breathers 
which are extremely short and steep envelope soli- 
tons, containing just one- two waves. 

Thus, though BFI is definitely a relevant parameter, 
but in some situations it does not completely charac- 
terize the freak waves, as it takes place, for example, 
for ID waves [H, [ll|, [H, [ll], or for ver y lo ng-crested 
waves [3, mi, or in crossing sea states [H, [l^. The 
reason is that coherent wave structures depend on addi- 
tional parameters as well. In the present work, we in- 
vestigate this question in more detail for weakly crossing 
sea states. More specifically, we consider sea states with 
two spectral peaks centered at wave vectors ko — Ak/2 
and ko + Ak/2 in the Fourier plane, and we assume 
|Ak| <C |ko|. Such a situation corresponds to the pres- 
ence of relatively long and wide wave stripes obliquely 
oriented to the wave fronts in a range of angles near the 
angle between the vectors Ak and ko- We study nu- 
merically how the process of rogue wave formation de- 
pends on the angle 6. The computations are based on 
the approximate theoretical model for long-crested fully 
nonlinear water waves, developed in Refs-jH, [2I] and 
later successfully applied in Refs.fsl, The model is 
intermediate between the exact Eulerian dynamics and 
the approximate equations for wave envelopes [general- 
izations of the nonlinear Shroedinger equation (NLSE)] 
suggested in Refs.[25|, [H, [S^. It should be emphasized 
that our method makes possible to compute profiles of 
individual waves, while in the works |l6l , only wave 
envelopes in crossing sea states for ^ = 7r/2 were studied. 

The main results of this work are the following. If 
< arctan( 1/^/2), then the nonlinearity is defocusing 
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along the stripes and it is focusing across them. The sit- 
uation is opposite at nearly perpendicular kg and Ak. 
Accordingly, rogue waves, occurring at the sea surface, 
are different in these two cases. In the first case, freak 
waves look as fragments of structures which are similar 
to the solitonic solutions of a focusing NLSE for a wave 
envelope. Such extremely narrow and steep solitons are 
in essence rows of alternating oblique crests and troughs 
(see Fig(T^). When perturbed by a weak two-dimensional 
(2D) random field, the oblique solitons can exist for many 
wave periods almost unchanged, but later they transform 
to zigzag structures similar to that described in Ref.[i] 
(see Figs|l]l2]for example). It should be noted that in the 
limit ^ ^ such extreme oblique solitons coincide with 
the recently discovered ID giant breathers [lO]. Thus, 
the fundamental role of these coherent structures in the 
dynamics of water waves is confirmed. At nearly perpen- 
dicular ko and Ak, another kind of coherent structures 
comes into play, similar to the dark solitons of a defo- 
cusing NLSE. Dark solitons develop at the interference 
minima and they transversally separate wave groups sub- 
jected to the longitudinal modulational instability. Freak 
waves in this case have nearly ID profiles, but they are 
bounded in the transversal direction by two dark solitons. 

To understand better numerical results, it is useful to 
have in mind a qualitative model describing weakly non- 
linear water waves in terms of a complex wave amplitude 
74(xi,X2,t), which determines the free surface elevation 
as follows, 

y'(xi, X2, t) ^ Re [A{xi,X2^t) ex.p{ikoXi — iuot)] , (2) 

where xi and X2 are horizontal coordinates, with xi along 
ko, ooq = {gk^Y^'^ is frequency of the carrier wave, and 
g is the gravity acceleration. The function A{xi,X2,t) is 
known to approximately obey a 2D NLSE [18], 



t=17 min 26.5 s. Surface elevation [m] 



i dA 
dt 



i dA 
2ko dxi 



(3) 

The oblique stripes roughly correspond to the following 
ID reductions of Eq.®: 

A = k^^^i^.r), ^ = ko[{xi-Vgrt) cos 0^X2 sin 0], (4) 

where r = c^o^, and Vgr = (cc;o/2A:o) is the group velocity. 
As a result, we have a ID NLSE describing the transver- 
sal dynamics of idealized, infinitely long wave stripes, 

i^r = \ [(l/2)cos2 6>-sin2 6>] + (5) 

Depending on the sign of the dispersion coefficient ol(S) = 
[(1/2) cos^ — sin^ 0], this is either focusing equation or 
defocusing one, and the dynamics is quite different in 
each case. For example, in the focusing case (when a > 
0), the nonlinearity can become saturated with the so- 
called (bright) solitons. 
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cosh[{s/V^){^-^o)] 



exp(-zrs^/4 + z^o), (6) 
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FIG. 1: (Color online). Evolution of a perturbed high- 
amplitude oblique soliton into a zigzag structure. 
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FIG. 2: (Color online). Maximum and minimum elevation of 
the oblique soliton vs. time. 



where s is a wave steepness, and ^o, (/>o are arbi- 
trary constants. These solutions describe infinitely long 
wave ridges consisting of alternating oblique crests and 
troughs. Physical conditions of applicability of the above 
formula imply 5 < 0.1 and s/^/a <^ 1, but actually these 
solutions have been found to continue qualitatively to 
considerably higher values s ^ 0.27. We specially stud- 
ied a long-time behaviour of such extreme solitons, both 
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FIG. 3: (Color online). Maximum elevation of the free surface 
vs. time in the numerical experiments A1-A4. 
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FIG. 4: (Color online). Experiment A3: the two big waves 
are at x ~ 1.6 km, q ^ [3.7 • • • 3.9] km, and at x ~ 1.5 km, 
[0.1 •••0.3] km. 

for 6 = (the giant breathers at 2D surface), and for 
0. It is one of the main results of the present work 
that in two dimensions extreme solitons can exist for a 
long time before transformation into zigzag structures. 
An example of evolution of a perturbed high-amplitude 
oblique soliton is presented in FigsHJO foi" ^0 ^ 100 m, 
^ arctan(l/5), and s ^ 0.22. 

In the defocusing case (when a < 0), the so-called dark 
solitons are possible, 

^ds = stanh [{s/^/^){^ - ^0)] exp(-zr5^/2 + i^o), 

(7) 

which separate two domains of opposite amplitude. 

In view of the above, it is clear that since the ef- 
fective dispersion coefficient a{0) changes the sign at 
0^ = arctan(l/V^), in the full 2D dynamics of random 
wave fields there should be two substantially different 
regimes, one regime at 6> < and another at close to 
7r/2. This hypothesis is confirmed in general by numeri- 
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FIG. 5: (Color online). Maximum elevation of the free surface 
in the numerical experiments B1-B3. 
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FIG. 6: (Color online). Experiment B2: the rogue wave is at 
X ^ 1.2 km, q ^ [1.0 • • • 1.4] km. 



cal experiments reported here. 

The computations were performed in the dimensionless 
square domain 27r x 27r with periodic boundary conditions 
along the horizontal coordinates x and q (see Refs.0,0 
for details). Thus, all the discrete Fourier modes corre- 
spond to integer wave vectors k = (/c,m). The vector 
ko was generally taken slightly different from the direc- 
tion of X axis, in order to take into account the effect 
of gradual re-orientation of wave crests along the oblique 
stripes (see Ref.[9]). Final results were rescaled to give a 
convenient for presentation value Aq ^ 100 m, which is 
quite typical in natural sea conditions. The correspond- 
ing wave period is Tq = [27tXo/ g]^^'^ ^ S s. Two small 
sets of typical numerical experiments are presented, des- 
ignated as A1-A4 and B1-B3. Within each set, at t = 
the normal Fourier modes of the wave field were taken in 
the form a/c^(0) = cF{k,m) e'xp{ijkm)^ with a positive 
function F{k,m) having two nearly Gaussian maxima at 
ko ± Ak/2, and with quasi-random initial phases jkm (in 
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FIG. 7: (Color online), (a) profiles of the freak wave from 
FiglU (b) 8 s later: "a hole in the sea" ; (c) 16 s later: the big 
wave has risen again. 



set A and in set B the phases ^ykm were not the same). In 
each experiment the coefficient c was different, thus re- 
sulting in different values of the total energy E. In set A 
we chose ko = (40.0, —2.5) and Ak = (7.0, 2.0), so a case 
< 0^ was simulated, while in set B it was a crossing sea 
state with = 7r/2: ko ± Ak/2 = (39.5, ±3.5). 

For set A, some results are presented in Figsl3]|4j The 



modulational instability acts in this case from the very 
beginning, and it needs a short time 5.. 8 min to produce 
freak waves in the initially most tall wave groups. The 
two neighbouring big waves in FigH] look as a fragment 
of an oblique soliton (compare to Figjl^). The computa- 
tions A4 and A3 were terminated at the moments when 
the freak waves broke, while in experiments A2 and Al 
the waves remained smooth, so at later times nearly sta- 
tionary, long oblique solitons were observed (not shown). 

Results of experiments B1-B3 (see Figsl5]|7j) more 
intriguing, since there were two stages in the evolution of 
the wave field before rogue waves arose. In the first stage, 
for 4.. 7 min after the beginning, it was the formation of 
dark solitons along interference minima, which process 
was accompanied by substantial decreasing of the wave 
amplitude along interference maxima. At the end of this 
stage, the free surface has been divided by dark solitons 
into domains of nearly ID dynamics. In the second stage, 
adjacent domains interact in a complicated manner, and 
in one of them the amplitude increases, resulting in fast 
development of the longitudinal modulational instabil- 
ity. As the result, a single rogue wave grows, which is 
squeezed from the lateral sides between two dark soli- 
tons, as shown in Fig|6l The rogue wave is "breathing" , 
with time-alternating tall crest and deep trough, and it 
approximately repeats the profile after 2To (see Fig I?]). 
After a dozen of the oscillations, the big wave spreads in 
the transversal direction and disappears (not shown). 

As results A and B are compared, it becomes clear that 
the unusual properties of the abnormal waves in weakly- 
crossing seas are more prominent when the interference 
stripes are nearly perpendicular to the wave crests. 
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